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1 Introduction 

Many physical systems display a strong dependence on the space dimension- 
ality, the best known example being given by phase transitions in equilibrium 
statistical physics. As for non-equilibrium systems, hydrodynamic turbulence 
shows a remarkable dependence on the space dimension as well. In three di- 
mensions, the kinetic energy flows from large to small scales in the form of 
a Kolmogorov-Richardson cascade. Conversely, in two dimensions, energy is 
transferred upscale at a constant rate, in an inverse cascade process [18]. Ad- 
ditionally, three-dimensional turbulence is characterized by a breakdown of 
scaling invariance and small-scale intermittency [13] , whereas the inverse cas- 
cade is apparently self-similar [5] and even shows some intriguing signatures 
of conformal invariance PQ . 

These observations have spurred the search of a critical dimension be- 
tween d = 2 and d = 3 in the hope that it could provide a starting point for 
an analytical attack of three- or two-dimensional turbulence, or both. This 
approach has been mainly applied to simplified models of turbulence, such 
as EDQNM [12], the shell model [TB], or a model obtained by generalizing 
the form of the two-dimensional stream function [24] . In those studies, the 
spatial dimension has been most conveniently reduced to a formal parameter 
that could take arbitrary values. The approach undertaken in the present 
paper differs from previous work at least in two important aspects. First, 
we shall consider a geometrical, rather than formal, way of looking in be- 
tween dimensions. Namely, we shall study the dynamics of fluid particles on 
cylindrical manifolds where the compactified dimension can be collapsed or 
inflated at will so as to connect continuously the two extreme cases. Second, 
we shall focus on a system that is fully under analytical control, that is the 
Kraichnan ensemble of velocities rather than Navier-Stokes turbulence. 

To study the turbulent transport of a passive scalar, Kraichnan intro- 
duced a Gaussian ensemble of decorrelated-in-time velocity fields [15] . A 
compressible generalization of the Kraichnan ensemble in the d-dimensional 
Euclidean space has been investigated by Gawedzki and Vergassola under 
the assumption of statistical isotropy [15] Fl In this model, the dynamics of 
fluid particles depends on three physical quantities: the space dimension, the 
degree of compressibility, and the (spatial) Holder exponent of the veloc- 
ity. The Holder exponent £/2 is greater than zero and less than one. This 
property mimics the behavior of a turbulent velocity field, whose realizations 
are typically non-Lipschitz in the limit of infinite Reynolds number. For any 
given d < 4 and < £ < 2, Gawedzki and Vergassola have identified a critical 
degree of compressibility separating two different phases of the Lagrangian 
dynamics. Below the critical value (incompressible or weakly compressible 
velocity fields), fluid particles separate superdiffusively. The probability dis- 
tribution of fluid-particle separations does not have a stationary limit in 
this regime. Above the threshold (strongly compressible fields), Lagrangian 
trajectories tend to collapse to zero distance, and the distribution of the sep- 
arations degenerates into a Dirac delta function. For d > 4, the former regime 



1 The smooth limit of this model had been previously considered in ref. [3]. 
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is the only possible one and the phase transition docs not occur|3 The above 
results have been subsequently elaborated by Le Jan and Raimond in the 
context of non-Lipschitz stochastic differential equations |21[[22] , 

Here we consider a generalization of the Kraichnan ensemble on a cylin- 
drical surface. A d-dimensional cylindrical surface can be constructed by 
taking R d and compactifying d — d! dimensions. The radius of the cylinder 
is the size of the compactified dimensions. When the radius tends to infinity 
we recover R d ; when it tends to zero we obtain R d . Thus, varying the ra- 
dius of the cylinder produces a smooth transformation from dimension d to 
dimension d' . 

We define a zero-mean Gaussian velocity field on a cylindrical surface by 
imposing the form of its covariance. We require that the covariance of the 
field tends to the one of the isotropic d-dimensional Kraichnan ensemble as 
the radius of the cylinder tends to infinity and to the one of the isotropic 
rf'-dimensional Kraichnan ensemble as the radius vanishes. The degree of 
compressibility is such that the velocity is weakly compressible in the limit of 
infinite radius and strongly compressible in the opposite limit. It is therefore 
possible to gradually move from one regime to the other by varying the radius 
of the cylinder. 

As we shall see, if in the limit of infinite radius the Holder exponent is 
equal to £, then in the limit of vanishing radius it is equal to £ + (d — d'). 
Hence, if attention is restricted to non-smooth velocities, the model under 
consideration is meaningful only when a single dimension is compactified 
(d! = d — 1). For the sake of simplicity, we shall conduct the analysis in two 
dimensions, where the two extreme cases are the two-dimensional plane and 
the straight line. We shall show that the dynamics of fluid particles results 
from two opposite effects. At small separations, Lagrangian trajectories ex- 
hibit a superdiffusive dynamics owing to the weakly compressible nature of 
the small-scale velocity. At large separations, fluid particles experience the 
trapping effect of a strongly compressible field. Consequently, the probability 
distribution of the two-particle separation tends to an invariant measure. This 
behavior is to be contrasted with the one observed in the two-dimensional 
isotropic case with the same Holder exponent and degree of compressibility. 

In the present context, the separation vector between two fluid particles is 
a stochastic process solving an Ito stochastic differential equation with non- 
Lipschitz diffusion coefficient. To guarantee the existence and the uniqueness 
in law of the solution, we shall add pure diffusion to the velocity field. By 
considering an appropriate Lyapunov function, we shall demonstrate that 
there exists an invariant measure for the fluid-particle separation. The in- 
variant measure is unique, ergodic, and non-degenerate as a consequence of 
the irreducibility and the strong Feller property of the process. 

The paper is divided as follows. Section[2]describes a generalization of the 
Kraichnan model on a d-dimensional cylindrical surface. The two-dimensional 
case is studied in detail in section [3l Sections |4] and [5] contain the results on 
the fluid-particle separation and its invariant measure. The limit of vanishing 
diffusivity, the effect of a viscous regularization of the velocity field, and the 
role of the Prandtl number are discussed in section [5] 



2 For d = 4, the collapsed phase can exist only for smooth velocity fields (£ = 2). 
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2 Kraichnan model on a d-dimensional cylindrical surface 

We consider the dynamics of fluid particles in a turbulent flow on a (i-dimcn- 
sional cylindrical surface S. The velocity field is a family of white noises 
taking their values in the space of vector fields on S. Specifically, v(t, x) is a 
Gaussian stochastic process with zero mean and covariance 

E(v a (t,x)vp(s,y)) = D af} (x - y)S(t - s), (1) 

where x,y G R d ' x [-ttL, TiL) d ^ d ' C R d {d > d') and L is the radius of 
the cylinder. The velocity is by definition statistically homogeneous in space, 
stationary in time, and invariant under time reversal. Moreover, we assume 
periodicity in the d — d! "radial" coordinates. 

It is convenient to write the spatial covariance D a p(r) in terms of its 
Fourier-space representation: 



with k = (k',k") G R d ' x ±1 d - d \ r = (r',r") G M d ' x [-ttL, nL) d - d ' , 
and a, (3 — 1, . . . , d. The presence of a series in the ^''-coordinates accounts 
for the periodicity of the velocity field in the r"-coordinates. 
We adopt the following form for the spectral tensor: 

F a0 (k) = ^ (fc;P) d+ , (2) 



with i G E+, ^ G [0, 2], p G [0, 1], and 

A a p(k; p) = (1 - p)5 aj3 + {pd - 1) 



As we shall see in the latter part of this section, F a p(k) has been chosen 
in such a way that, in the limits L — > and L — > oo, D a p(r) tends to the 
covariance of an isotropic random field|f| 

The spectral tensor is real, symmetric (F a p(k) = Fp a (k)) and non- 
negative definite Vfc G R d ' x \ J J d - d ' , i.e., 

F af} (k)u a u > V(u u ...,u d ) G K d , 

l<a,/3<d 

as can be checked using the Cauchy-Schwartz inequality. These properties 
guarantee that D a p(r) is the spatial covariance of a homogeneous random 
field (e.g., ref. [25], p. 20). Moreover, F a p(k) is an even function of k, and 
therefore the velocity is statistically invariant under parity: D a p(—r) = 
D a0 (r). 

3 The spectral tensor could in principle be multiplied by a positive coefficient 
determining the intensity of the velocity fluctuations. For the sake of simplicity, we 
set that coefficient to one. 
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As a consequence of statistical homogeneity and parity invariance, the 
covariance of velocity differences can be expressed in terms of D a p(r): 

E([v a (t,x + r)- v a (t, x)][vp(s,x + r) -v p (s,x)}) = 2d aj3 (r)5{t - s) (3) 

with d a p{r) = D a0 (O) - D aP (r) [25]. 

The meaning of the parameters p, £, and £ may be understood by con- 
sidering the limit of D a p{r) for L — > oo and for L — > 0. 

The limit L —> oo (and l/L-> dfc") yields: 



Urn D^(r) = — ^ / dfc" f dk 

L ^°° {2w) a J R d-a> J R d> (\\k'\\ 2 + \\k"\\ 2 +£- 2 ) 



= * / dfc e^(fc; f ) 
(2?r) 7 R d (||fc|| 2 + ^- 2 )~ a ~ 

In this limit, D a p{r) tends to the spatial covariance of a d-dimcnsional 
isotropic field with correlation length £ and degree of compressibility p [15] 
ITT] . The parameter £/2 represents the inertial-range Holder exponent of the 
velocity: Y? a=1 d aa {r) = 0{\\r\\t) as \\r / £\\ -► 0. For £ = the velocity field 
is purely diffusive; for £ = 2 it is spatially smooth, and its spatial regularity 
decreases with decreasing £. In particular, the Kolmogorov scaling is obtained 
for £ = 4/3, for the time integral of must be proportional to ||t|| 4 ^ 3 in 
Kolmogorov's phenomenology [13J0 

It is worth noting that for a finite L equation (HJ) describes the velocity 
covariance at space separations much smaller than L. 

In the second limit, L — > 0, we obtain^ 



zfc' -r' 

(||fc'|| 2 +^- 2 )^ 



with £' = £ + (d - d') and 



iC=l + ^l ifrf>l, 



p(d-d') + d'-r d'-l 
A X x{k';p') = 1, £=p(d-l) if d' = 1 and p > 0. 

4 The same conclusion can be reached rigorously by defining the Kraichnan en- 
semb le a s the limit of an Ornstein-Uhlenbeck process for vanishing correlation 
time [TO] . 

5 This can be shown by multiplying D a p(r) by a function f(r"), integrating 
over r" £ [— tvL,ttL), taking the limit L — *• 0, and noting that only the term 
corresponding to k" = has a non-zero limit equal to 



IT r ik' r' 
/(0) 7^7 / , dk' * 



A a0 (k';p') 



(llfe'll 2 +l~ 2 ) i ^ L - 
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We thus recover the covariance of a d'-dimensional isotropic velocity field 
with Holder exponent £'/2, correlation length £, and degree of compressibil- 
ity p' [M]. 

The exponents £ and £' must satisfy the inequalities < £ < 2 and 
< £' < 2. Therefore, the limit L — > makes sense only in two cases: 

a) d' = d-l, £e [0,1], and £' = £+1 € [1,2]; 

b) d! = d - 2, £ = 0, and £' = 2. 

We are interested in the situation where the fluid particles disperse when L — > 
oo (cZ-dimensional isotropic flow) and collapse when L — > ((/'-dimensional 
isotropic flow). Moreover, we focus on spatially rough velocity fields leaving 
aside the cases £ = and = 2. This situation can be realized only in 
case a), for £ € (0, 1), and under the conditions [T5] : 

d , . d 1 

p < — and p > -j. 

The first inequality is actually satisfied for all d and £ given that £ € (0, 1) 
and p G [0, 1] . The second inequality can be rewritten in terms of p as follows: 

d-2 

p -WTYy (5) 

The restriction < p < 1 and inequality {5| imply the additional condi- 
tion d < 5. 

In the remainder of the paper, we shall investigate the statistics of fluid- 
particle separations on a two-dimensional cylindrical surface (d = 2). 



3 Two-dimensional cylindrical surface 

For d — 2, case a) is the only realizable one, corresponding to d' = 1. Con- 
dition (0 reduces to p > independently of £. 

The spatial covariance of the velocity field takes the form 

D afj {r) = —j- > e t / dfci 6 

47T^L ' J R r, 2 , (o\ 2 , 1 1 

3 = -oo [k ± + [j) + F J 

with r = (ri, r2) E f2 = Ix [— 7rL, 7rL). In eq. ((6]) we have written k = (k%, fe) 
with ^2 = i/L and j € Z to make the dependence on L explicit. We shall 
keep this notation in the remainder of the paper. 

We now restrict attention to space separations much smaller than I. For- 
mally, this is equivalent to considering the limit £ — > oo. The spatial variance 
of the velocity field, D a p(0), diverges as I tends to infinity (appendix P£f : 
this behavior reflects the divergence of the average kinetic energy of the 
fluid. Nevertheless, d a p(r) has a finite limit for all r, and the statistics of 
velocity differences remains well defined. 
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The limit of d a p(r) for £ — > oo can be computed explicitly (appendix [A}. 
The correlation of the axial component is written: 



lim dn(r) = 



p 








2 3 +«7T 


3/2/- 


(i + 4)i 



In I 



+ 



# g r i- e ( l + (l-pX r / r i + ; 



27r 3 / 2 r (2 + 4 



2 J J'=l 



2 cos 



J r 2 

L 



P 1 + 



2L 



+(l-2p) 



2£ 



#3 



(¥) 
,(¥) 



(7) 



where K v (z) denotes the modified Bessel function of the second kind of or- 
der v and argument z. The correlation of the radial component has the form: 



lim d 2 2(r) = 



(1-p) 


r 




i 


2 3 +«7r 3 / 2 r 







In I 



1+5 



2^(2 + |) 1 



_ £ j i + pj r f i + € 



2 cos 



J r 2 

L 



(1 - p) 1 + 



+(2p-l) 



n /0>ij 



2L 



2L 



(8) 



Finally, the mixed correlations can be written as follows: 



lim di 2 (r) = lim d 2 i(r) 



2^ r (, + i )U' WWW 



(9) 



The limit ^ — > oo will be hereafter understood. 
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3.1 Large-scale form of the covariance of velocity differences 

To understand the nature of the random velocity field, it is useful to consider 
the covariance of velocity differences at space separations much greater than 
the radius of the cylinder. 

The series in eqs. (JT])-© converge uniformly (appendix [A"f. For |ri|/Z — > 
oo, it is therefore possible to deduce the asymptotic expansion of d a( a(r) 
from the limiting behavior of the single terms of the series. The asymptotic 
expansion of K v {z) for z — > oo is (e.g., ref. [SJ, formula II 7.13(7)) 

Kv{z) ~{r z e ~ z (|argz| <37r / 2 )- 

Thus, the ^-dependent contributions to d a p(r) decay exponentially fast with 
increasing space separation. The remaining contributions give 



with 



and 



d u (r)~SiN 1+s + xi as^^oo (10) 



Si = 



p 


r 


;- 








2 3+5 7r 3/2 r 






L 



*i = r \^ L C(i + 0- 



In the latter equation 

CM = E? ( s > !) 

is the Riemann Zeta function. Likewise, we have 



d 2 2{r) ~£> 2 |ri| 1+ « + k 2 as ^ -► oo (11) 

with 



p 1 + (1 - p)i 

Finally, the off-diagonal terms vanish at large space separations: 

lim d 12 (r)= lim d 21 (r) = 0. (12) 

\ri\/L — >oo \r\\/L — >oo 

The above asymptotic expressions show that, at separations much greater 
than the radius of the cylinder, the velocity difference may be regarded as 
the superposition of two independent one-dimensional random fields. One 
field is directed along the axial direction; the other one is directed along the 
radial direction. Both the fields depend only on r%. In particular, the axial 
field is a one-dimensional Kraichnan velocity field with Holder exponent 1 + £ 
like the one considered in ref. [3D]. 

At large separations, the small-scale dynamics manifests itself through an 
effective diffusivity represented by the constants xi and >t 2 . 
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4 Fluid-particle dynamics 



In the present context, the separation between two fluid particles can be 
regarded as a stochastic process on Q with diffusion coefficient d a p(r) (and 
drift coefficient equal to zero) . To ensure the (weak) existence and uniqueness 
of the trajectories of the process, we add diffusion to the velocity field and 
replace d a p{r) by 



n>0. 



The additional term can model the action of molecular diffusion on fluid 
particles as, e.g., in ref. [15j . The constant k will be referred to as diffusivity. 

The separation vector between two fluid particles will be denoted by R. 
According to the above remark, R satisfies the Ito stochastic differential 
equation^ 

dR(t) =V2a(R(t))dB(t), R{0)=ref2, (13) 

where B is Brownian motion on Q and a is defined through the Cholesky 
decomposition of the matrix d K : 



d K 



with 



a(r) 



fVd^) 
dUr) 




d? x (r) J 



r e Q. 



\VdUr) 

Although the diffusion matrix d K is nor Lipschitz continuous nor bounded, 
the existence and uniqueness of the solution of eq. (I13| can be proved using 
Stroock's and Varadhan's theory of martingale problems [28]. To directly 
exploit this theory, we shall first consider the periodic extension of eq. (fT3|) 
on R 2 , and then project the resulting process on fl. We therefore introduce 
the projection p : M 2 — > fl with 



P(r) = ( 



T\ , —nL + ri — 2ttL 



T2 

2ttL. 



era . 



and define a := a o p. Likewise we denote d K := d K o p 

Before proceeding further, it is convenient to define some notation. The 
spaces of bounded measurable and bounded continuous functions on fl will 
be denoted by £§b{£i) and respectively. The set ^ 2 (f2) will be the 

space of functions having two continuous derivatives. Analogous definitions 
will apply to functions defined on M 2 . 



6 If X and y denote the positions of two fluid particles, the separation vector 
between the two particles is defined as R := Y — X . The common physical notation 
for the evolution equation for R would be 

where £ is white noise and the statistics of 5jtv(t) := v(t,Y(t)) — v(t,X(t)) is 
defined by eq. © 
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Proposition 1 The ltd stochastic differential equation on R 2 : 

dR(t) =V2a(R(t))dB(t), R(0) = reR 2 , (14) 

where B is the standard Brownian motion on M 2 , has a unique (in law) weak 
solution. In particular, the solution is a continuous Markov process. 
For reR 2 and UCl 2 measurable, let 

P(0, r; s, U) := P (r(s) € U if R(0) = rj 

be the transition probability distribution of R, and let {Tt) t>0 be the associate 
transition semigroup: 

Ttf(r) := [ f(p)P(0,r;t,dp) 

with f £ ^b(R 2 ). Then, the semigroup {Tt) t>0 has the strong Feller property, 
i.e. f t (M b (R 2 )) C %(M. 2 ) for all t > 0. 

Proof The diffusion matrix has the following properties: 

1. d K is continuous; 

2. d K (r) is symmetric and strictly positive definite for all r 6 R 2 . The dif- 
fusivity k is indeed assumed to be strictly positive, and the spatial co- 
variance of velocity differences must be symmetric and uniformly non- 
negative definite for all r (ref. [25] , P- 97), i.e., 

d a p{r)u a Uf3 > V (iti, it 2 ) and r S R 2 ; 

l<a,/3<2 

3. there exists a positive constant C\ such that for all a, /3, and r 

\fr a p{r)\<C x {\ + \\vf). (15) 

This property is a consequence of the asymptotic behaviors ((TO]) to (fT2]) 
and of the the fact that d K {r) is bounded at the origin and continuous 
on R 2 . 

Under the above conditions, Stroock's and Varadhan's uniqueness theorem 
apply to the martingale problem for d K 28J. Then, the proposition follows 
from the equivalence between the well-posedness of martingale problems and 
the existence and uniqueness in law of weak solutions of stochastic differential 
equations (ref. [27], pp. 159 and 170). □ 

The process R can be regarded as the projection of R on Q: R(t) = p(i?(t)) . 
The properties of R can then be deduced from those of R. 

Corollary 1 Equation (|13p has a unique (in law) weak solution. The tran- 
sition semigroup (T t )t>o associated with R has the strong Feller property: 
T t {m h (n)) c %{f2) for all t > 0. 
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Proof Given an initial condition r € fl, a weak solution of eq. (113j) can be 
constructed by taking a solution of eq. (fT4l) with the same initial condition 
and projecting it on fl. 

The key observation to prove uniqueness in law is that any solution of 
eq. (|13|) on fl can be uniquely mapped into a continuous solution of eq. (|14p 
on R 2 . Then, uniqueness in law in R 2 guarantees that also the solution on fl 
is unique in law. 

Finally, (T t ) (>0 has the strong Feller property and the projection p is 
locally invertible (with continuous inverse). Hence, (Ti) t >o has the strong 
Feller property. □ 



5 Invariant measure of fluid-particle separations 

An invariant measure for R can be constructed by adapting to the case 
under examination the procedure described in ref. |26] . Clearly, an invariant 
measure may exist only if the trajectories of the stochastic process do not 
"escape to infinity" . To control the behavior of the first component of R, 
which is not bounded, we therefore introduce the Lyapunov function V : 
fl -> R+: 



V(r) = 



h{ h + l)c 4 + 2(1 - h 2 )c 2 r 2 - h(l - h)r\ 

4(1-/0^+1) it|n|<c, 
1 inl 2 ^) if H>c, 



2(l-/i) 



where c > and 1 > h > 0. The function V is twice continuously differcn- 
tiable and has the asymptotic behavior needed for the proof. 

Lemma 1 If (1 +£)/2 > h > 1/2 and A denotes the infinitesimal generator 
of (T t ) t>Q , then the Lyapunov function has the following properties for all £ € 

(o,i).-" 

1. limjj^i^ AV(r) = -oo; 

2. there exists m G R such that AV(r) < m for all r £ fl; 



3. T t V{r) = V(r) + / T s AV(r)ds. 
Jo 

Proof For / S ^(fl), the infinitesimal generator of (Tt) t>0 has the form 
Af(r)=tiia(r)a T (r)D 2 f(r)}, 

where D 2 f denotes the Hessian of the function /. The action of the genera- 
tor A on V(r) is written: 

From eq. (fT0|) we obtain 

AV{r) - -(2ft - l)D 1 |ri| 1+s - 2ft as ||r|| ^ oo. 



12 



If (1 + C)/ 2 > ft- > 1/2, wehavelim|| r ||^ 00 .AF(r) = -oo for all £_ £ (0, 1). It is 
worth noting that this latter result relies on the fact that dn(r) = 0(\ri\ 1+ ^) 
as ||r|| — > oo with 1 + £ > 1. 

Property 2 is a consequence of the continuity of AV and of property 1. 

Finally, the transition semigroup satisfies 



T t f(r) = f(r)+ [ T s Af{r)ds 
Jo 



for any / G ^?(/2). The same property holds true for the function V, as can 
be shown using an approximation procedure similar to the one described in 
ref. pp. 167-168. The details are given in appendix iBl □ 

We now make use of the properties of the Lyapunov function to obtain the 
following result. 

Proposition 2 There exists an invariant measure fi for the stochastic pro- 
cess R. i.e. 



n 



(T t f)(r) f,(dr) = /_/(r)/i(dr) 
for all f £ and for all t > 0. 

Proof Given j-q £ f2, we have 



(16) 



T s (m -AV) (r )ds = m 



m + 



T s AV(r )ds 
V(r ) - T t V(r Q ) 



<m + 



V(r ) 



Hence 



sup 

t>i 



T s {m - AV){r )ds 



< oo. 



(17) 



We now introduce the family of "average" measures (nt)t>i on Q defined, 
for any / € %(.!?), as 



1 



f(r)fM(dr) = - I TJ(r )ds. 



We show that ([it)t>i is uniformly tight. For a given N > we define 

E N = {re Q:m- AV(r) < N}. 

The set En is compact: it is closed since it is the preimage of a closed 
subset of Q, and must be bounded since AV(r) — > — oo as ||r|| — > oo. As a 
consequence of Markov's inequality, the measure of the complement of En 
satisfies: 

ME C N ) < ^J [m-AV{r)]fH(dr) = jf t J Q T s (m - AV)(r )ds. 
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Using eq. (fl"7]) we conclude that for all e > there exists En C fl with 



N = - sup 

e t>i 



T fl (m-^)(r )d4 



such that fit(E N ) < e for all i > 1. The family (fit)t>i is therefore uni- 
formly tight. Then, there exists a measure /i and a sequence (t p ) p >o with 
limp^oo t p — oo such that /i tp converges weakly to \i as p — > oo. This means 
that J n fd[j,t p — ► Jq f d[i as p — * oo for all / € (e.g., ref. [Bj, theo- 

rem 11.5.4, p. 404). 

We now show that fi is invariant. For any / £ %(J3), for t > 0, and for 
all p such that t p > t, we have: 



/(r-) Mip (dr)- / T t f(r)fx tp (dr) 



-/ T s f(r )ds~- T s T t f(r )ds 



p Jo 



1 



P ^0 

1 



T s f(r ) ds 



T s f(r )ds 



pt-\-tp 



T s f(r )ds 



p Jo 



t+t p 



TJ(r ) ds 



< 



2*||/|| 



Hence, for alH > 0, 

f(r)(H p (dr) 



T tf{r) IM p (dr) -> as p -> oo. 



The semigroup (T t ) t >o satisfies T t (^b(/2)) C since (T) t >o has the 

strong Feller property. By using the weak convergence of fit to p, we can 
thus conclude that ||TBJ) holds for all / E %(fi) and for all t > 0. The 
measure is therefore invariant for R. □ 

To show that the invariant measure is actually unique, we need the following 
result stating that R has no closed invariant set different from the whole 
space. 

Lemma 2 The semigroup (T t )t>o is irreducible, i.e. the transition probabil- 
ities of R, P(0,r;t,U), are strictly positive for all t>0, for all r £ fl, and 
for all non-empty open sets UC.fl. 

Proof For all r G M 2 the linear transformation associated with cf(r) is in- 
vertible, and therefore maps M 2 into itself. Hence, the semigroup (Tt) t>0 is 
irreducible (ref. [29], theorem 24, p. 66). 

The transition probabilities of R arc connected to those of R as fol- 
lows: P(0,r;t,U) = P{0,r*;t,U*) where U* = p^{U) and r* is any point 
in p _1 ({r}). Therefore, (T t ) t > is irreducible. □ 



We can now state the main result regarding the invariant measure of R. 
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Theorem 1 There exists a unique invariant measure /i for the stochastic 
process R. The measure fi is ergodic and equivalent to any transition proba- 
bility P(0,r;t,U) with r £ Q, t > 0, and U C Q measurable. Moreover, fi is 
absolutely continuous with respect to the Lebesgue measure, and is therefore 
non degenerate (i.e., broad in r). 

Proof We have already proved that fi is invariant. Its uniqueness, ergodicity, 
and equivalence to any transition probability follow from the fact that the 
transition semigroup associated with R has the strong Feller property and is 
irreducible [7]fT7] (see also ref. [5], chapter 4). 

To prove the absolute continuity of [i with respect to the Lebesgue mea- 
sure, we introduce the family of transition probabilities 



/•OO 

Q A (r, U) = A / e- As P(0, r; s, U) ds 
Jo 



with r £ Q and U C fl measurable, as well as the associate transition 
semigroup 

Txf(r)= f f{y)Qx{r,dy). 
J a 

Likewise, we define an analogous family Q\(r,U) for the process R. The 
measure fi is invariant also for (7\)\>o: 

f (i(dy)T x f(y) = \ f dse- Xs f ^(dy)T s f(y) = f n(dy)f(y) 
J n Jo Jo Jn 

for any / £ 38b(fi), and hence 

lx{U)= f v(dy)Qx(y,U) (18) 

for any measurable set U C Q. 

For all r £ R 2 , the measure Q\(r,-) is absolutely continuous with re- 
spect to the Lebesgue measure (see ref. [25], theorem 10, p. 24). It fol- 
lows that Q\(r,-) has the same property for all r £ Q since Q\(r,U) = 
Q\(r* ,p^ 1 (U)) for a given r* £ p~ 1 {r}. From eq. (fT5|) , fx is therefore abso- 
lutely continuous with respect to the Lebesgue measure. □ 

As a consequence of lemma [2J and theorem [T] the transition probability of R 
has a positive density with respect to the Lebesgue measure: P(0, r; t, dp) = 
p(0, r;t, p)dp. The probability density function is the (possibly weak) solu- 
tion of: 

d t p = Mp, (19) 

where, for / £ tf 2 (Q), 

Mf(p)= J2 d Pa d p ,d^(p)f(p). 

l<«,/3<2 
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6 Conclusions 

We have studied the dynamics of fluid particles in a compressible turbulent 
velocity field on a cylinder. The model that we have introduced is a general- 
ization of the isotropic Kraichnan ensemble. Although the parameters of the 
velocity have been set in such a way as to produce explosive separation of the 
fluid particles in the isotropic limit (L — > oo), on the cylinder the probability 
distribution of the separation tends to an invariant measure. This behav- 
ior is a result of the compressibility effects generated at large scales by the 
compactification of the "radial" dimension. 

The diffusivity k has been taken strictly positive to guarantee the exis- 
tence of solutions to eq. (|13|) . The addition of Brownian motion to Lagrangian 
trajectories influences the dynamics of fluid particles at small separations. 
Therefore, the presence of a nonzero diffusivity may be relevant for the non- 
degeneracy of the invariant measure, but should not affect the existence of 
the invariant measure itself, which rather depends on the large-scale form 
of the velocity field. The limit n — > may be tackled by means of Wiener 
chaos decomposition methods 20,21,22,(53]. We conjecture that our results 
remain valid in that limit. Indeed, in the situation considered, the small-scale 
dynamics of fluid particles is the same as in the weakly compressible phase 
of the isotropic Kraichnan ensemble. In that regime, Lagrangian trajectories 
separate in time even for vanishing k owing to the poor spatial regularity of 
the velocity |15j . Thus, the invariant measure should remain non-degenerate 
as k — ► 0. 

The Reynolds number is infinite in our study since the viscosity of the 
fluid, v, has been set to zero from the beginning. For the same reason the 
Prandtl number Pr = v/k is equal to zero. The viscosity can be taken 

. . 2 II jr ||2 

into account by multiplying the spectral tensor (J2J) by the factor e~ n 11 11 , 
where r\ cx v 3 ^ 4 plays the role of the viscous length of the flow [TT|. This 
modification has a small-scale regularizing effect on the velocity field, which 
for any rj > is locally Lipschitz continuous. Obviously, a positive r\ does not 
alter the proofs of the results shown in the paper. 

The order of the limits k — ► and rj — > 0, however, deserves a detailed 
discussion. Taking the limit k — > before 77 — ^ is equivalent to letting Pr 
tend to infinity. The opposite order corresponds to the limit Pr — ► 0. As 
first observed in ref. [SJ, when these limits are considered the range of weak 
compressibility splits into two ranges: what is now called the range of weak 
compressibility in the strict sense, < p < (d — 2 + £)/(2£), and the range of 
intermediate compressibility, (d— 2+£)/(2£) < p < d/^ 2 . In the former range, 
the order of the limits k — > and rj — > is not relevant for the Lagrangian 
dynamics [14] . At small scales, fluid particles disperse irrespective of the order 
of the limits, and therefore we expect the invariant measure of the separation 
to be non-degenerate. By contrast, the order matters in the latter range [14] • 
For intermediate values of the compressibility, if the viscous regularization is 
removed before the diffusivity (Pr — > 0), the small-scale Lagrangian dynamics 
is once more characterized by the explosive separation of the trajectories. If k 
goes to zero before rj (Pr — > 00), the trajectories coalesce also at small scales, 
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and the invariant measure of the separation should degenerate into a Dirac 
delta function. 

In summary, we believe that the present study captures the behavior of 
the Lagrangian trajectories on cylindrical manifolds for all Pr except for the 
limit Pr — > oo in the intermediate-compressibility regime. These results are, 
moreover, relevant to turbulent transport of passive scalar fields in virtue of 
the relation subsisting between the scalar correlations and the dynamics of 
fluid particles [MET]. 

We conclude by noting that when both the dimensions of the plane are 
compactified one obtains the Kraichnan flow on a two-dimensional periodic 
square studied in ref. 4j. The velocity field considered there was however 
smooth in space. 

Acknowledgements The authors are grateful to F. Flandoli and Y. Le Jan for 
fruitful discussions. 



A Covariance of velocity differences 

For d = 2, d' = 1, £ (E (0, 1), f £ (1, 2), the spatial covariance of the velocity field 
takes the form 



i r 



L- 

with r £ Q. We first establish the convergence of the above series. It is convenient 
to denote the integrals by D^(n) and thus rewrite the covariance as follows: 

i=— oo 3=1 

(20) 

Using the inequality 

\A a p(k); p)\ < 1 - p+ \2p- 1 Va,/3 = 1,2 and Vfe £lx yZ, 

Li 

we obtain that the coefficients of the series satisfy for all t\ 

P$(n)| < (l-p+|2 P -l|)Af, 
with (e.g., ref. [9], formula I 1.5(2)) 



h \ l 2 p) r i + 4 Vi 2 ^ 2 



The series Mi converges for all £ £ (0, 1) and £ > as well as in the limit I — » 

oo, as it can be checked by means of the integral test. Th en, the Weierstrass criterion 
guarantees that the series in the right-hand-side of eq. (|20[) converge uniformly and 
absolutely on £2. The uniform convergence will allow us to compute lim^oo d a p{r) 
by exchanging limit and summation. 
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The basic analytical ingredient to derive eqs. J?])-© is 

° lfeiri 2 1 - J -K 1/2 K v {\zr 1 \) \ri\" 



(21) 



with Re{v) > -1/2 and | argz| < vr/2 (e.g., ref. [5], formula II 7.12(27)). 

We first compute the correlation of the axial component of the velocity; the 
correlation of the other components may be easily derived from Du(r). 

In the limit I — > oo, we have 



lim Dli /0) (ri) 



= (1 - 2P) ^ 



ifcl7- 1 



dki 



dki 



[(l-p)£ + p*?] 



i- 2 -i- it 
K i + 7^ 



+ p dki 



kfe ikiri 



(fc? + £) 



2+i 



ifciri 



+ p / dki 



ifeiri 



(l-2p) 



+ P 



7r l/2| ri |(3 +e )/2 K ^ ( |i r , i | ) 

2(i+«/ 2 r (2 + 1) 

rfiTT) ' 



To compute D^' (0) we can use the asymptotic expansion of K v (z) for z — > 



^ (!)" - 



2 

r(v) (x 



2 



2) 



Hence we obtain 



lim JJ<f O) (0) = p 



20,-1) 

j -^^ 2 r(i±i) 



(1/ < 1) (22) 
(|)"' y + 0(^) (Kv<2). (23) 



r 1 



+ (1 - 2p) 



-1-5 



^ 2 r( 3±i) 



r 2 + - 



-i-£ ^/ 2 



^ 1/2 ( i + g-pQr(i±£) 
2r (2 + 1) 



For j = and £ < 00, we have 



0?(n) = p 



o ifc l r l 



dki 



2~vr 1/2 .. .1+1 T , nri\\ 



{kl+l- 



The asymptotic expansion (|22p shows that 13^ (0) diverges like as £ — > 00. By 
using expansion (|22[) . it is nonetheless possible to show that 



lim [D^(0)-D^(n)]=p 



(0), 



^ 1/2 |^(-^)1 

2 1 +«r(i + |) 



ri 
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(note that r( — (1 + £)/2) < for < £ < 1). Hence, the covariance of the axial 
component of the velocity difference, dn(r), has a finite limit as £ — > oo. 
For the other components we have: 



^(n) = 



dfcl 



[(l-p)fc? + p^] 



(*? + £) (*? 



3 -4- 



2+C 
1 \ 2 



and 



/o 



dfci 



fcie ifeiri 



Therefore, ^^'(ri) can be derived from (n) by replacing p with 1 — p: 

2( l +o/ 2 r (2 + 1) 
i+i 2 ( 1 -«/ 2 7r 1 / 2 J ft: i±i (|f n|) |n|( 1+ «)/ 2 



hm D(f°)(r 1 ) = (2p-1) 



+ (1 - P) 



lim£><f o) (0) = 



-i-t 



hm [l?g ) (0)-i?W(ri)] = (l-p) 



(0)/ 



2 r( 2 + |) ' 

ff Va |r (- i±«) 



2i+cr(i + §) 



The mixed correlation can be obtained, for j 7^ 0, by differentiating formula (|21[) 
with respect to ri and by using -$-\x v K v {x)\ = — x v K u -\(x) (e.g., ref 9 , for- 
mula II 7.11(21)): 



L L 



2^ i r(2 + |) 



*I±£ (li»-i|)sgn(r-i). 



Hence 



lim D^ 0) (0) = 0. 



L»f 2 ) (r 1 ) = Vn€ 



For j — we have 

Finally, eqs. (JTJ)— (|9j| may be derived by recalling that 

1 00 



and using the uniform convergence of the series in the right-hand-side of eq. (I20|l . 
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B Proof of Lemma Q] 

To prove property 3 of lemma [1] we first observe that 

E {\\R(t)\\ 2 ) = \\R(0)\\ 2 +2E (J \\a(R(s))fds} < C 2 (l + J E (||i?(s)|| 2 ) d^j , 

where Ci > and ||cr|| := [Tr(acr T )] 1 ^ 2 . The above inequality is a consequence 
of (|15[) . Gronwall's inequality then yields 

E(||ii(t)|| 2 ) < C 2 c C2t (24) 

for all t > 0. 

Following ref. [26], we consider, for all 7 > 0, the function ip~, : R+ — > K+ with 

Jz 0<z< 7 
^ (2) ~ Wt+1) 7 + 1 < 2 

and ip 7 G < ^' 2 (R+) and monotonically non-decreasi ng. Moreover, we define V 1 := 
tpy o 1/. Applying .4 to V~, and taking into account (|10[) yield 

|AVr(r)| = c&(r) 



with 3 + £ - 4ft < 2 as ||r|| — > 00. 

Since V 7 G %f 6 2 (i?), we have for all 7 > 

TtV^r) = Vy(r) + [ T s (_4V 7 )(r)ds. 
Jo 

We now show that each term of the above equation tends as 7 — » 00 to the corre- 
sponding term in property 3 of lemma [T] 

Obviously, ]inL y ->< x Vy(r) = V(r) for all r G Q. Likewise, Vy(R(t)) / V{R(t)) 
almost everywhere as 7 —> 00. Therefore, by the monotone convergence theorem 
limbec T t Vj(r) = T t V(r). Finally, lim 7 ^cx> AVy{R(t)) = AV(R(t)) and, from 
eqs. Q25) and ®, |AVy(JlC*))| < C 3 (l + ||-R(i)|| 2 ) with C* 3 > and E (\\R(t)f) < 
00. Then, it follows from the bounded convergence theorem that 

lim f T a {AV^)(r)ds = [ T s {AV)(r)ds. 

T^°°io JO 

This concludes the proof. 



( dV 



O 



3+£-47i 



(25) 
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